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Abstract 



Exploring the idea of Broadhurst on the "sixth root of unity" we present an ansatz 
0> ' for construction of a basis of transcendental numbers for the e-expansion of single 
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Qh! three-loop master integrals are calculated. 
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1 Introduction 

One of the interesting cases of diagrams occurring in QCD and the Standard Model is of the 
propagator type with one mass and the external momentum on the mass shell (the bubble 
(vacuum) type diagrams also belong to this class). Only a limited number of analytical 
results is available in this case: the two- and three-loop diagrams occurring in QCD are 
given in Refs.[|l], |^; the analytical result for the two-loop bubble integral are presented up 
to the finite part in [Q, |^, Q and the 0{e) term in 0; the finite part of three-loop bubble 
master integrals is calculated in 0; the results for all other two-loop on-shell propagator 
type integrals with one mass are collected in |^ . For some diagrams the results are available 
in terms of hypergeometric functions H, 0. However, the e-expansion of hypergeometric 
functions is an independent and sometimes very complicated task. 

At the same time a lot of different methods for numerical calculations of Feynman dia- 
grams are available now |10|, |Tl|, [T2|, ^, |T^. From another point of view the development 
of numerical algorithms, in particular PSLQ [|1^, allows to determine with high confidence, 
whether or not particular numerical results for master-integrals can be expressed as a ra- 
tional linear combination of several given transcendental constants. However, the problem 
arises of how to define the full set of the transcendental structure (the "basis") occurring 
in given diagrams. The structure of massless multiloop integrals is well understood at the 
present time. In contains mainly the ^-function [|16| and Euler sums [l^. Recently the 



connection between knot theory and the transcendental structure of Feynman diagrams was 
established |jT^. In more complicated case, when dimensionless parameters exist, harmonic 
sums |T^ , pp] ] appeared. For studying the massive case, the differential equation method ^1\ 



and geometrical approach |Q Q can be useful. 

The aim of our investigation is to find the transcendental structure of single mass scale 
propagator type diagrams. In this paper we restrict ourselves to the consideration of di- 
agrams with two particle massive cuts only. All our results are obtained empirically by 
carefully compiling and examining a huge data base on high precision (about 200 decimals) 
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2 Basis 

Our aim is to construct a basis (the full set of transcendental numbers) for single mass scale 
propagator type diagrams with two-particle massive cuts. The simplest diagrams of this 
type are shown in Fig.|l]. 

An algebra for such a basis can be constructed as follows: the product of two basis 
elements having weights Ji and J2, respectively, must belong to the basis of weight Ji + J2. 
The weight of basis elements can be connected, e.g., with powers of e within the Laurent 
expansion of dimensionally regularized [Q Feynman diagrams. 

In constructing the basis of weight J we will use Broadhurst's observation 0, that the 
sixth root of unity plays an important role in the calculation of the diagrams (see also 
^^)- Accordingly we conjecture, that the basis of weight J contains functions obtained 



by splitting of polylogarithms Lij (z) , Lij (1 — z) into real and imaginary parts, where z 
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Figure 1: Examples of the simplest single mass scale integrals with two-particle massive 
cut. Bold and thin lines correspond to massive and massless propagators, respectively. 
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) can be reduced to these two 



exp u| j , exp [i^-] ■ The general case, z = exp (±«f + 

values |2^. E.g., with Li2„ f e'^j = Ghn {&) + ^Cl2n {9), the basis of weight 2n must contain 

Gbn (f) , Gl2n (^j , Cl2n (f j and Cl2n f^j- Not all of these functions are independent: 

Gl2„ (cutt) ~ C,2n-, where uj is an arbitrary rational number; (1 + 2-*^^^") Cl2n (^j = Cl2n (f )• 
The basis we are going to construct is supposed to hold for any number of loops and 
arbitrary order of the ^-expansion. Consequently we connect the lowest weight basis elements 
with Lii [z) and Lii (1 — z), which results in tt and In 3. Comparing this with the one-loop 
diagram of Fig.|l| (ONSll), the expression given in p^ allows us to find from its finite part 
the transcendental normalization factor 4=, so that our lowest basis is given by ^ and In 3. 
In general the normalization factor 4= is introduced according to our experience. 

Using the above conjectures, we obtain the basis of higher weights: for weight 2, e.g., we 
have -%ln3,C2, ~m^^2 (f) and ln^3. The results up to weight 5 are presented in Table I. 

Some transcendental structures like, e.g., Lss (| j do not appear in the table since they could 
be eliminated by means of relations given in p6[. Including weight 3 the relations found in 



26[| were sufficient. For weight 4 we needed the new relation 
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and for weight 5 we used the further relations 
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in order to restrict the basis as much as possible. 
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The above table contains the following functions, which appear in the splitting of the 
polylogarithms into real and imaginary parts [^ . 
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The basis of weight 5 is sum of columns 5'* and 5^ in Table I. Some numerical values of these 
functions for special values of their argument are given in appendix A. 

It should be noted, that the number Nj of basis elements with weight J satisfies the 
following relation: 



Nj = 2N 



j-i 



2\ 



(4) 



3 Examples 

Let us present several examples of two- and three-loop diagrams expressible in terms of 
the basis structure elaborated in the previous section. Here we will use the notation 152 = 
-4^Cl2 (?) and are working in in Euclidean space-time with dimension A^ = 4—26. Moreover, 



the normalization factor 



''(in)^ 



r(i+e) 



for each loop is assumed. 



As first example we consider the two-loop single mass scale bubble diagram (VLlll in 
Fig.|l]). Taking into account the 'magic connection' ||^ between the one-loop off-shell massless 



triangle and two-loop bubble integral and using the one-fold integral representation p^ for 
the former, we have: 
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(5) 



Another form of the integral representation of this diagram can be found in ||^ . In this way, 
the calculation of the coefficients of the e-expansion is reduced to the following integrals: 



1 



dx 



X + x"^ 



ln"xln''(l -x + x"^ 



(6) 



The values of these integrals for a + b < 2 are know from [0, ^, ^ . The results for the next 
two indices (a -|- 6 = 3, 4) are given in appendix B. Combining all these results we have: 
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This expression coincides with analytical results obtained before. The e -part was obtained 
in [^ and all higher terms were obtained by Davydychev p8| . 

Let us consider another example: JOll diagram of Fig.|I| with the external momentum on 
the mass shell of the internal mass. The first two terms of the e-expansion of this diagram has 
been calculated analytically in []^. Here we present the next two terms of this expansion. It 
is well known that two diagrams with different sets of indices are needed as master integrals 
For the numerical investigation we choose the integral with indices 122 (the first index 



belonging to the massless line) and a linear combination investigated in ||30|. To obtain the 
numerical results for these two integrals with high precision |^ we have constructed their 



Fade approximants from the small momentum Taylor expansion [13|. The Taylor coefficients 
were obtained by means of the package [^. Within the above mentioned basis we have: 



and 
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Figure 2: Three-loop master-integrals. Bold and thin lines correspond to massive and 
massless propagators, respectively. 
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The values of the master integrals (with indices 111 and 112) can be obtained by means of 
recurrence relations [16|| implemented in the package ON-SHELL2 [33|. They are relatively 
lengthy and therefore will not be presented here. 

Now we consider master integrals occurring in the package |^: 1^3(0, 1,0, 1, 1, 1, 
£'3(1, 1, 1, 1, 1, m) and .05(1, 1, 1, 1, 1, 0, m) as shown in Fig.^ the notation being as in 



m 



These integrals contribute in calculations like [p5[| . 

Results for D3(0, 1,0, 1, 1, 1, m) and i?3(l, 1, 1, 1, l,m) have not been published before. 
Two independent approaches have been applied to obtain high precision numerical values 
for these integrals. One method started from analytical results |jl|, 0, ^ for two-loop prop- 
agator diagrams. The finite three-loop diagrams were then obtained by multiplication with 
the massive propagator with power j > 1, resulting in a one-fold integral representation, 
which was integrated by MAPLE with accuracy more than 200 decimals. Alternatively the 
£^3-integral was calculated by the method which was used earlier for the numerical calcu- 
lation of 1)3(1, 1, 1, 1, 1, 1) [^. Here the product of the two one-loop propagator diagrams 
were expanded in their small/large external momentum and the integration was performed 
piecewise. Details of such calculation were published in P7[ 
decimals was also obtained in this manner. 

The analytical results are the following ones: 



An accuracy of about 100 
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Let us now consider the diagram ^5(1, 1, 1, 1, 1, 0) shown in Fig.^. As shown by Broad- 
hurst PI it contains the new transcendental number Vs^i. We were able to verify this, but 
we also found the following decomposition in terms of our basis elements: 

^l-.^^! ,.V2^^ , 13,_,^, 259, , 3^ n)/27r 



^3 1 = —S' - -ttLss — + — Cs In 3 C4 + -Ls^'^ — . (12) 

^' 16 ^ 4 V 3 / 24^'^ 108^ 8 ^ V 3 / ^ ' 

4 Conclusion 

In this work we have constructed the basis of transcendental numbers for single mass scale 
diagrams with two-particle massive cuts. 

The basis , up to weight 5, is given in Table I. The conjectures which have been used 
for the construction of the basis are as follows: the set {hj} of transcendental numbers of 
weight J has the structure: 

hj = {hj-KhK, 6j}, ir = 1, 2, ■ ■ ■ , J - 1, 
where new elements {6j} are connected with functions arising in the splitting of polyloga- 
rithms Lij (2;), Lij (1 — z) into real and imaginary parts, and z = exp u| j , exp [i^j ■ Up 
to the weight 5 we observed also, that the the number Nj of basis elements with weight J 
satisfies the relation (H). 

The fact that the 'sixth root of unity' [^] plays such a crucial role in the construction 
of our basis remains mysterious. The introduction of the 'normalization factor' 4= remains 
a matter of experience. Nevertheless our construction allows us to calculate several new 
diagrams (or obtain new terms in the e-expansion) in an analytical form. 

The interesting question remaining beyond our work is how to construct the basis for 
single-mass scale diagrams having several different massive cuts. From multiloop QED 
and QCD calculations P, |^, |l^, |3^, ^ we know that other transcendental numbers, like 



TT, In 2, G, Li„ (I) or their products appear. In this case the structure of the basis elements 

is connected with polylogarithms of 2; = exp u|j ,exp{m) @]. At the present moment we 
don't know how the mixing basis could look like. 
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APPENDIX A. Numerical values of basic elements. 

Here we present numerical values up to 40 decimals of functions, which for special values 
of their argument provide our basic elements. These numerical values were obtained by 
numerical integration of their integral representation by means of MAPLE. 
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Ls3 f ^) = -2.1447672125694943914826421122153512974711 ■ ■ ■ , 

-- 6.0094975498188889162047887062032707405970 ■ ■ ■ , 
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-- -0.5181087868296801173472656387316967550219 ■ • • . 



APPENDIX B. Analytical results of integrals occur- 
ring in the e-expansion of the two- 
loop bubble. 

Here we present the values of the integrals (|^). Their numerical values were obtained by 
MAPLE with accuracy of 230 decimal and again the analytical results were obtained by 
PSLQ. 
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